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Abstract - The finite momentum transfer (q) longitudinal optical response a L (q,co) of graphene 
has a peak at an energy uj — TiVFq. This corresponds directly to a quasiparticle peak in the spectral 
density at momentum relative to the Fermi momentum ftF — q. Inclusion of coupling to a phonon 
mode at cue results, for w < |we|, in a constant electron-phonon renormalization of the bare bands 
by a mass enhancement factor (1 + A) and this is followed by a phonon kink at uie where additional 
broadening begins. Here we study the corresponding changes in the optical quasiparticle peaks 
which we find to continue to directly track the renormalized quasiparticle energies until q is large 
enough that the optical transitions begin to sample the phonon kink region of the dispersion curves 
where linearity in momentum is lost in the renormalized Dirac Fermion dispersion curves and the 
correspondence to a single quasiparticle energy is lost. Nevertheless there remains in a L (q,uj) 
features analogous to the phonon kinks of the dispersion curves which are observable through 
variation of q and u. 
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Important information on the charge dynamics of the 
Dirac fcrmions in graphene is obtained in optical absorp- 
tion experiments. Results are reviewed by Orlita and 
Potcmski. [T] For bare bands the real part of the zero 
momentum limit optical conductivity, a(q = 0, uS), has 
a Drude peak around lu = which has its origin in the 
intraband optical transitions. There is also an additional 
piece due to the interband transitions with onset at twice 
the chemical potential which provides a constant univer- 
sal background of value a = ire 2 /2h. [2H3 In experi- 
ments [5] the Pauli blocked region between the Drude and 
the onset of the universal background is not found to be 
small as expected in a bare band model but instead has 
a value of order <7o/3. This observation can be partially 
explained due to impurities, electron-phonon inter- 

actions and/or coulomb correlations, although its precise 
origin remains controversial. In conventional metals the 
electron-phonon interaction is known to renormalize many 
of their properties. [TU] It leads to incoherent phonon- 
assisted Holstein side bands [TT] in addition to the main 
coherent Drude response with optical spectral weight now 
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distributed between these two parts. The effect on the 
Drude weight can be expressed in terms of the electron 
phonon mass renormalization parameter, A. The Drude 
weight is reduced with a multiplicative factor of 1/(1 + A) 
while the remaining A/ (1+A) is transferred to the sideband 
due to the new absorption processes in which a phonon is 
created by a photon along with a hole-particle pair. In 
graphene, electron-phonon effects have been seen [H] in 
the electronic density of states (DOS), N(u), measured in 
scanning tunnelling spectroscopy (STS). These observa- 
tions are expected in systems where the DOS is energy 
dependent. P^HTB] Phonon 'kinks' have also been seen 
in angle-resolved photoemission spectroscopy (ARPES) 
[T7JHH] as predicted theoretically. [TI5] Recently, near-field 
optical techniques have been used to obtain information 
on the finite momentum conductivity |20H23| rather than 
the standard long wavelength, q — > 0, limit. This has al- 
lowed for nano-imaging [2"0H2"2"] of the graphene plasmons 
and could in principle be employed to get information on 
plasmarons [23] a scattering resonance of an electron and 
a plasmon. |24j In this paper we show how a{q,uj) can 
provide information on electron-phonon renormalization 
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Fig. 1: (Color online) The real part of the finite momentum 
longitudinal optical conductivity, a L (q,u))/o~o as a function of 
oj/ho for a selection of values of q/ftF from 0.0 — > 0.9 (as 
defined in legend). These include an electron-phonon inter- 
action with A ~ 0.18 for a single optical phonon mode at 
oje = 200 meV (shown by vertical dashed line, where here 
Ho = 1 eV, and include a residual scattering r\ = 0.005/io. In- 
set: — ImE EPI (u))/ no as in Eq. @ for parameters in the main 
frame. 



effects complementary to ARPES and STM. 

The Kubo formula for the real part of the finite momen- 
tum optical conductivity is given in the simplest bubble 
approximation, which neglects vertex corrections but in- 
cludes interactions through a self energy, by 



a(q,uj) 8 

CT oj 
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where A s (k,uj) is the spectral function of the Dirac 
fermions. The scattering overlap factor, F ss i ((/>), has the 
form F ss > ((p) = i [1+ss' cos </>] where is the sum of the an- 
gles of initial momentum, k, and final momentum, k + q, 
with respect to the k x axis (zigzag direction). Also, we 
define the angle a for q and obtain the relevant overlap 
factor 



FsA4>) = \ 



, fccos(26») + qcos(6 + a) 

s . 

^Jk 2 + q 2 + 2kq cos(6» - a) 



(2) 



In this representation, 9 is the angle over which k is inte- 
grated in Eq. (TTJ). The values of s and s' are each either 
+1 or —1 which represents the upper (+) and lower (-) 
Dirac cones and for simplicity wc will take the chemical 
potential to fall in the upper cone, representing an electron 
doped sheet. In general the conductivity, <Jij(q,oj), can 
be written as a linear combination of a longitudinal part, 
a L (q,oj), and a transverse part, a T (q,oj). If we select 



q = qx, where a = in our orientation, then one would 
find that the conductivity along the x-direction, a xx (q 1 oj) 1 
has only a contribution from the longitudinal conductiv- 
ity (ie. a xx (q,oj) = a L (q,oj)) while the conductivity along 

the y-£LXlS, ®yy 

(q,oj), has only a contribution from the 
transverse part (ie. a yy (q,oj) = a T (q,oj)). Throughout 
this work, we are interested in examining quasiparticle- 
likc peaks which appear only in the longitudinal part of 
the near-field optical spectra, and therefore will limit our 
discussion to examining a L (q,oj). 

Denoting the self energy by T, s (k,oj) the spectral func- 
tion is 

1 |ImE s (fe,a;)| 



A s (k,oj) 



n[uj- Re£ s (fc,o;) - e|] 2 + [ImS s (fc,w) 



2 ' 



(3) 

where ef. = svpk — fx with vf the Fermi velocity. For 
the electron-phonon interaction (EPI) the self energy docs 
not depend explicitly on momentum or band indices but 
changes with frequency. For coupling to a single phonon 
at frequency u>e, the imaginary part of the self energy, 
Y, EPI (ui), is given by 

— ImS EP/ (a), cj.e) = 



w — We + no\, for uje < oj < W c 



IT A 

— \0J + UJE + Ho 



for 



UJE > OJ > 



- Ho + OJE, 



Ho - OJE j 



(4) 

from which the real part can be obtained through a 
Kramcrs-Kronig transformation. We use the notation set 
out in Ref. [5] where more details can be found. In Eq. (fj[|) , 
A is a constant that can be adjusted to get a desired value 
of the mass enhancement parameter, A, and W c is a cut- 
off on the bare band energies adjusted to get the correct 
number of states in the Dirac approximation of two val- 
leys in the Brillouin zone. By definition A is obtained 
from the small oj limit of the real part of E(w). In this 
limit, ReT, EPI (oj) = ReS £P/ (0) - Xoj where the constant 
piece, Re£ BP/ (0), shifts the chemical potential from its 
bare to interacting value. Detailed calculations of the 
electron phonon interaction in graphene have been pre- 
sented by Parks and co-workers [2oTI?7] who conclude that 
a reasonable approximation to the complete calculations 
is to use a model of coupling to a single phonon of en- 
ergy oje = 200 meV and here we follow this suggestion. If 
one were inclined, a distribution of phonon modes can be 
created through an integration over oje in Eq. H]weighted 
by the phonon density of states at each frequency, as has 
been done previously. [51l2"8"] 

The real part of the longitudinal conductivity follows 
by a numerical evaluation of Eq. ( [TJ . Results are shown 
in Fig. [1] for ten values of q/kF from 0.0 — s- 0.9. The 
solid black curve is the well known 6 9 q = case and is 
included here for comparison. There is a Drude piece at 
small oj which is followed by a boson assisted part that sets 
in abruptly at oj = cue- Here we have used a model for the 
self energy which also includes a constant impurity scatter- 
ing term, 77, in addition to the self energy scattering term 
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Fig. 2: (Color online) Schematic of doped Dirac bands. 
An electron-phonon interaction renormalizes the Dirac bands, 
changing the slope through the Dirac point and producing 
kinks at uj — ±uje where here, the Fermi level (w = 0) is 
illustrated by a change in color. 



from the electron-phonon interaction. This constant scat- 
tering contributes no energy renormalization but broadens 
the Drude into a Lorentzian. The EPI however, reduces 
the optical spectral weight of the coherent Drude part and 
redistributes it in the phonon assisted part which corre- 
spond to intraband absorption. At higher energies (not 
shown here) there is additional absorption due to inter- 
band optical transitions which in the bare band case start 
abruptly at a frequency of twice the chemical potential. 
With electron-phonon coupling these also spill into the 
Pauli blocked region below 2/i. When q is finite we get 
the series of curves labelled by the energy of the position 
of the peaks in a L (q,uj). One might refer to these peaks 
as quasiparticle peaks since their energy would correspond 
exactly to e^ F _ g in the bare band case. 

Here for small values of q, the peaks in <j(q,u>) cor- 
respond instead to the renormalized quasiparticle ener- 
gies e + (kp — q)/[l + A]. However, when q approaches 
the phonon energy this simple renormalization ceases to 
hold. The broadening of the curves increases abruptly at 
w = we- This abrupt onset indicates that a new channel 
for scattering has opened; in our case, phonon emission. 
In the inset of Fig. [TJ we plot the imaginary part of the 
self energy, which contains a sharp jump at the phonon 
energy. The broadening due to phonons is clearly seen 
upon comparing results for q < 0.2 and q > 0.2. As q is 
increased the broadening increases as expected from the 
inset of -ImS(w). Note that the Dirac point in the inset 
occurs where the imaginary part of the EPI self energy is 
zero, which comes from a zero in the electronic density of 
states at this energy. 
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Fig. 3: (Color online) The peaks in a L (q,u)) plotted as a func- 
tion of kp — q coincide with peaks in the EPI renormalized 
spectral function, A(k,ui). For comparison, the bare disper- 
sion is shown by the dashed line, while the green line with 
circles is the bare dispersion divided by (1 + A) which is the 
slope of the EPI renormalized spectral function (shown in red) 
at k = kp. The peaks in a L (q,uj) sample the renormalized 
dispersion until near uj = — uje- Inset: same as main frame but 
showing both positive and negative frequencies. 



Fig. [2] illustrates the effect of the electron-phonon in- 
teraction on the Dirac cone dispersion relations. For the 
bare case, the Fermi velocity sets the angular dimension 
of each cone. Interactions distort the cones in two impor- 
tant ways. First of all, around the Fermi level the bare 
Fermi velocity is renormalized to a dressed value through 
a mass enhancement factor of (1 + A). As we move away 
from the Fermi energy this simple law starts to break down 
and phonon structures (referred to as 'kinks') develop at 
uj = ±uje from the Fermi energy, and are measured in 
ARPES experiments. These features are shown in the in- 
set of Fig. [3l The solid red curve depicts the peaks in the 
phonon renormalized spectral function which illustrates 
the kinks at ±cj as compared to the dashed black curve 
which is the bare cone dispersion. In general the renor- 
malized dispersion curves are also broadened. 

Here we want to know how the phonon 'kinks' in the 
dressed dispersions manifest in the near field optics. This 
is illustrated in the main frame of Fig. [3] where the ener- 
gies of the peaks in the real part of the optical conductivity 
a(q, uj) are traced as a function of momentum, kp — q, and 
are shown as the solid black curve. These are to be com- 
pared with the solid red curve based on the peaks in the 
spectral density, A(k, uj) as a function of momentum, rela- 
tive to the Fermi level. The open circles with green curve 
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Fig. 4: (Color online) The real part of the finite momentum 
optical conductivity, (J L (q, oj) /ao as a function of uj/ho for a se- 
lection of values of q/kF- The dashed lines include the electron- 
phonon interaction identically to Fig. [T] for q/kF displayed in 
the legend, while the solid lines include only the residual scat- 
tering of rj — 0. 005/jo but for shifted q values (marked near the 
peaks). The solid circles mark the analytical clean limit re- 
sults, with an adjusted Fermi velocity to obtain the same peak 
positions. 



represent a straight line having a Fermi velocity renor- 
malized by 1/(1 + A) as compared to the bare dispersions 
shown as the dashed curve. Note that for energies a bit 
above the phonon energy at lue the peaks in the optical 
conductivity and in the single particle spectral density of 
ARPES agree with each other and both fall on the renor- 
malized Dirac fermion energy e£/(l + A) (light solid green 
curve). Above cue this simple relationship ceases to hold 
with optics (solid black) roughly tracing the renormalized 
dispersion (green curve) while ARPES (solid red) is closer 
to the bare dispersion (dashed). While the dressed spec- 
tral peaks show a phonon kink at w = lue there is no 
discernible structure in the optical conductivity at this 
energy. It is important to remember that optics involves 
a joint density of states from an initial occupied electron 
state to a final unoccupied one. While many optical tran- 
sitions correspond to a given energy it is clear that tran- 
sitions from a kink to a second kink have a large weight 
and thus introduce structures into a L (q,ui) for frequencies 
which fall in magnitude just below |2u;b|. 

In Fig. [3] we concentrated on the energies of the renor- 
malized peaks in the optical conductivity and their rela- 
tionship to peaks in the spectral density. Both are due to 
the real part of the quasiparticle self energies. In Fig. @] we 
illustrate instead the effect of interactions on the broaden- 
ing of the peaks which depends instead on the imaginary 
part of the self energy. For each of the four values of q/kp 
we show 3 cases. The dashed lines include both real and 
imaginary parts for the full self energy including coupling 
to the phonon as well as to impurities (?y) . Open circles are 
for the clean limit (no scattering) case but where the ener- 



gies e^T of the bare bands have been shifted to peak at the 
same energy as the renormalized bands (values labelled 
on individual peaks). The solid curves are similar in this 
regard, but include the residual scattering, r\ = 0.005/^o- 
Starting with the large q case, peaked near uj — 0.82/io, we 
see the broadening near the peak in both the residual scat- 
tering and phonon cases in comparison to the clean limit 
which drops to zero above to = q. Thus there is substan- 
tial filling of the Pauli blocked region between lu = q and 
lu = 2fi — q due only to scattering. Note the difference be- 
tween residual scattering and phonon cases. In this region 
there is an order of magnitude difference in a L {q, w). This 
is in contrast to the same curves as lu — > 0, where the two 
cases are nearly identical, while the clean limit again drops 
to zero. In the low frequency limit — ImS £P/ (w) « 0, leav- 
ing only the value in the residual case. Thus in this limit, 
separation between the dashed (residual) and solid (EPI) 
curves must be due to energy renormalization effects which 
come from the (1 + A) factors of the real part of the EPI 
self-energy. By contrast for the case of q < uje, illustrated 
in the blue curves labelled at lu = 0.08/io, all three deviate 
much less from each other in the region of the peak, and 
again, it is not until frequencies above lu = lue that large 
distinctions begin to occur. 

Very recently, Ashby and Carbotte [29] have studied 
the case of residual scattering alone and for small q val- 
ues have derived in first order an analytic formula which 
has proved quite accurate. Here we have generalized their 
work to include the electron-phonon interaction in the 
regime where the renormalized quasiparticle energies are 
well represented by the formula e,t/(l + A). We obtain 



a (q,u) _ 4/xrj 

<7 7T 



2cos 2 (60(2?y) 



2tt (2t?) 2 + [lu(1 + X) + q cos ( 



(5) 



In the clean limit (/; — » 0) we get 

cr L (q, lu) 8/jq ( ^ 
n V q 



ctq 



yjq 1 — cj 2 1 + A ' 



(6) 



with q = q/(l + A). We see that the square root singularity 
is at lu = q rather than at q as in the bare band case. Also 
the optical spectral weight under a L (q, lu) is reduced by a 
factor of 1/(1 + A) and is given by 



doj 



a L (q,cu) 



(TO 



2/x /(l + A). 



We see that in the case of a finite value of rj, taking out 
the (1 + A) factor next to lu in Eq. (JSJ will leave an overall 
multiplicative factor of 1/(1 + A) in addition to changing q 
to q. The scattering rate is also renormalized in the same 
manner as fj = ?//(l + A). This result means that the 
electron-phonon interaction reduces the effective residual 
scattering rate at low frequencies. These scalings are all 
in agreement with previous work [5] for the conventional 
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case of q = which exhibits a Drude peak (intraband 
transitions) in the energy region below the inset of the 
interband transitions which start at w = 2(Iq- 

To conclude, we have calculated the effect of electron- 
phonon coupling on the finite momentum optical response 
of graphene with the aim of getting a first understand- 
ing of how this interaction modifies the bare band re- 
sults. A Kubo formula is employed in the simplest bub- 
ble approximation which neglects vertex corrections, but 
includes electron-phonon rcnormalizations in the Dirac 
fcrmion spectral density A{k,u>). In addition, there ap- 
pears in the formula for the finite q conductivity an im- 
portant scattering factor, F ss > ((/>), which incorporates chi- 
ral properties of the Dirac charge carriers associated with 
initial and final electronic states. As such, this factor de- 
pends in a fundamental way on the direction of q relative 
to the crystal lattice and direction of measurement. The 
conductivity in the ^-direction, a xx , can be decomposed 
for a general direction of momentum q into a linear combi- 
nation of the longitudinal response for q scattering which 
occurs along the zigzag direction (k x ) and a transverse 
part for q scattering along the armchair direction (k y ). 
Here, a peak is identified in a L (q,ui) which is directly re- 
lated to the quasiparticle energy £ + (/cf — q) in the bare 
band case. The electron-phonon interaction is found to 
shift the position of this peak in energy and to broaden 
it. Nevertheless, for values of q small enough that only 
optical transitions involving both initial and final states 
with energies well below the phonon energy uje enter, the 
peaks in a L (q, uj) still track perfectly those in the electron 
spectral function A(k, uj) measured in ARPES and both 
relate directly to the renormalized quasiparticle energies 
e + (kF — + A). However, as the magnitude of q is in- 
creased, the quasiparticle dispersion curve starts to devi- 
ate significantly from linearity and exhibit a phonon kink. 
At this point quasiparticle and optics begin to deviate 
from each other. For example, optics shows no kink struc- 
ture at u> — uje- Instead the phonon structure in a L (q,ui) 
is found to have shifted to higher energies in the region of 
uj < 2uie- In this energy range, optical transitions involv- 
ing a phonon kink in both initial and final states become 
possible and this leads to deviations of the optical spec- 
trum from linearity in analogy to the phonon kinks at uje 
found in the renormalized Dirac fermion energies. In this 
case there is a contribution from other optical transitions 
which emphasize less phonon structure and so the image 
of phonon is not as sharp in optics as it is in ARPES. 
Finally, we found that at higher energies, uj > uje, the 
optical peaks continue to track well the renormalized en- 
ergies e + (kE — <z)/(l + A) while the quasiparticle peaks 
move instead towards their bare band values. 
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